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Abstract 

We use first-principle calculations for transition metal impurities V, Cr, Mn, Fe, Co and Ni in 
GaAs as well as Cr and Mn in GaN, GaP and GaSb to identify the basic features of the electronic 
structure of these systems. The microscopic details of the hole state such as the symmetry and 
the orbital character, as well as the nature of the coupling between the hole and the transition 
metal impurity are determined. This could help in the construction of model Hamiltonians to 
obtain a description of various properties beyond what current first-principle methods are capable 
of. We find that the introduction of a transition metal impurity in III-V semiconductor introduces 
a pair of levels with t2 symmetry - one localized primarily on the transition metal atom referred 
to as Crystal-Field Resonance (CFR) and the other localized primarily on the neighboring anions 
referred to as the Dangling Bond Hybrid (DBH). In addition, a set of nonbonding states with e 
symmetry, localized on the transition metal atom are also introduced. Each of the levels is also 
spin-split. Considering Mn in the host crystal series GaN GaP — > GaAs GaSb, we find 
that while in GaN the hole resides in the t*^^^ level deep in the band gap, in GaAs and GaSb 
it resides in the t^^^ level located just above the valence-band maximum. Thus, a DBH-CFR 
level anticrossing exists along this host-crystal series. A similar anticrossing occurs for a fixed host 
crystal {e.g. GaAs) and changing the 3d impurity along the 3d series: V in GaAs represents a 
DBH-bclow-CFR limit, whereas Mn corresponds to the DBH-above-CFR case. Consequently, the 
identity of the hole carrying orbital changes. The symmetry (e vs. t2), the character (DBH vs. 
CFR) as well as the occupancy of the gap level, determine the magnetic ground state favored by 
the transition metal impurity. LDA-I-U calculations are used to model the effect of lowering the 
energy of the Mn d^ state by varying U. We find that this makes the DBH state more host-like, and 
at the same time, diminishes ferromagnetism. While the spin-splitting of the host valence band 
in the presence of the impurity has been used to estimate the exchange coupling between the hole 
and the transition metal impurity, we show how using this would result in a gross underestimation 
of the coupling. 

PACS numbers: PACS number: 75.50.Pp,75.30.Et,71.15.Mb 
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I. INTRODUCTION: ELECTRONIC STRUCTURE AS A TOOL FOR CON- 
STRAINING MODEL HAMILTONIANS 



The prospect of manipulating the electron spin to store and transport information in 
semiconductor devices has led to renewed interest in the physics of transition metal (TM) 



impurities in semiconductors - an area w 
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ferromagnetism (FM) at ambient temperatures has led to the investigation of the mechanism 
that stabilizes FM in transition-metal doped semiconductors. One useful approach to obtain 
an understanding of the electronic properties of these systems is the first principles electronic 
structure approach where one focuses on the explicit electronic and spin wavefunctions of 
the system. Variational minimization of the total energy, determines within the underlying 
approximations of the spin density functional theory some of the basic features of the states 
involved, such as the extent of localization, the magnitude of the spin interactions as well 
as the identity of disorder and compensating defects (antisites; interstitials) . However, the 
approach does have the drawback of underestimating the extent of electron correlations in 
addition to being a zero-temperature approach. While comparison with experiment {e.g. 
ferromagnetic temperature vs. alloy composition [13]) can be used to assess the extent to 
which electron correlations are underestimated, the first-principle results have to be generally 
mapped onto a model Hamiltonian to calculate finite temperature properties. 



While model Hamiltonians 
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the Hamiltonian are rarely justified in their own right. Generally [14, 
renormalizes away the electronic degree of freedom and retains only the spin degree of 
freedom for the transition metal impurity. One then assumes a local interaction between 
the transition metal impurity and the free carrier (usually RKKY-like), and then solves for 
various physical properties of these systems. 

In the present work we use first-principles calculations to examine whether the assump- 
tions made in model Hamiltonian treatments are consistent with an ab-initio description of 
the electronic structure of these systems. Our detailed results are then cast in the language 
of a simple electronic structure model, which could be used in an informed construction of 
model Hamiltonians. 
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We start by identifying the main physical quantities that come into play in determining 
the electronic structure of these systems. When a trivalent cation site such as Ga of a Ill-V 
semiconductor is replaced by divalent Mn, an acceptor level (denoted as E(0/-)) is generally 
created in the band gap. If the Fermi level ep lies below this E(0/-) level, then Mn is charge 
neutral, i.e., its formal oxidation +3 equals that of the Ga atom being replaced. In this 
case there is a hole in the Mn-related orbital. If, on the other hand, the Fermi level is above 
E(0/-), then the Mn-related orbital captures an electron from the Fermi sea {i.e. creating 
a hole there), becoming negatively charged {i.e. oxidation state Mn^+). In this case the 
hole resides in the Fermi sea. The Mn-induced hole for ep < E(0/-) features prominantly in 
contemporary theories of ferromagnetism 

The model Hamiltonians involve three entities - the host crystal states, the transition 
metal atom, and the impurity- induced hole state. There are approximations made at various 
levels which involve decoupling various degrees of freedom. At the first level, one decouples 
the orbital degrees of freedom associated with the transition metal atom, describing it with 
a localized spin-only part. The spin is interacting with a hole system through a local 
exchange interaction. At the next level of approximation, one reduces the problem to that 
of the transition metal spin interacting with the hole spin, assuming that the host crystal is 
unperturbed. The main assumptions made in such approaches, which we wish to examine, 
are: 

(i) The hole resides in a hulk-like, hydrogenic, delocalized state. This picture is based on 
the assumption that the perturbing potential Vnin{r)-VGa{r) generated by the impurity is 
dominated by a long-ranged Coulomb part, as a result of which only a small percentage of the 
charge resides in the Wigner-Seitz cell and the rest is distributed over a large portion of the 
host crystal. In this "host-like hole" picture, one reduces the problem to a quasi- hydrogenic 
form in which the acceptor state is designated via quantum numbers {s, p ...) of the host 
lattice. In such cases the wavefunction of the acceptor level is delocalized, and can essentially 
be constructed from the host crystal F states. This picture is motivated by the fact that 
divalent, post transition metal atom elements such as Zn^+ form in III-V semiconductors 
quasi-hydrogenic acceptor states [2^ with small binding energies. Similarly, extrinsic p- 
type doping of II- VI dilute magnetic semiconductor CdMnTe |2l| also form hydrogenic hole 
states. However, unlike Zn in GaAs or extrinsic p-type CdMnTe, the Mn atom introduced 
into III-V's has chemically active d orbitals so it is not obvious that the acceptor state 
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it forms in GaN, GaP or GaAs would qualify as a delocalized host-like hydrogenic state. 
Indeed, the microscopic features determining the localization of the hole wavefunction, such 
as the (i-character of the acceptor level must be considered. Such interactions could change 
the symmetry {t2 vs. e) of the hole state, hence its coupling to the host. The pertinent 
quantum designation of the hole state is impurity-like (^2, e) not host-like effective-mass 

(ii) The host valence band maximum (VBM) levels are unperturbed by the transition metal 
impurity. In this view, the host band structure represented in the model Hamiltonian could 
be described by a k.p model, valid for the pure host crystal and is decoupled from the part 
of the Hamiltonian involving the host -|- hole system. However, since one of the symmetry 
representations of the Mn d orbitals in tetrahedral sites (^2, e) is the same as that of the 
zincblende VBM (t2), such states could couple, hence become mutually perturbed. 

(iii) The spin of the hole couples to the spin of the transition metal impurity via an inter- 
atomic local exchange interaction Jpd- As only the spin degree of freedom of the transition 
metal atom is considered, while the orbital degree of freedom is ignored, the free carriers feel 
the effective magnetic field produced by the transition metal impurity spin. This is modelled 
as a local exchange interaction, Jp^, between the transition metal impurity and the spin of 
the free carrier. Hence, the magnitude of Jpd determines the energy scale of ferromagnetic 
ordering. Areas visited by the free carrier are rendered ferromagnetic. However certain ma- 
terials are found to show activated behavior in their transport implying no free charge 
carriers, yet they exhibit ferromagnetism. The current model which requires delocalized 
carriers cannot explain ferromagnetism in such systems. 

In what follows, we use first principles calculations to examine the validity of assumptions 
(i)-(iii) for "id impurities V, Cr, Mn, Fe, Co and Ni in GaAs, as well as for Mn and Cr in 
GaN, GaP and GaSb. We then construct a qualitative model that explains our numerical 
results. We find that 

(i) The Mn- induced hole could have significant 3rf character. The assumption of a "de- 
localized hydrogenic hole" is not supported by first-principle calculations. The depth of the 
acceptor level (reflecting its localization) and the coupling of the ?)d impurity orbitals to the 
hole change markedly with the host crystal in the series GaN GaP — > GaAs GaSb. The 
hole generated by introducing Mn in GaN is found to have significant 3rf character, while in 
GaSb the hole is found to have primarily host character. Further, the symmetry of the hole 
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depends on the combination of the host crystal with the impurity atom. For example while 
in GaAs:Co the hole has dominantly ^2 character, the corresponding isoelectronic impurity 
ZnSe:Fe has a hole with e symmetry. 

(ii) The presence of the transition metal impurity perturbs the valence band of the host 
semiconductor. The extent of the perturbation depends on the relative position of the 
impurity generated levels (referenced to the valence band maximum of the host) which have 
the same symmetry as the valence band maximum. In GaAs:V which has levels with e 
symmetry in the bandgap, the perturbation is small, while for GaAs:Mn the perturbation 
is large. 

(iii) The interaction between the spin of the TM atom and the spin of the host-like hole 
has a predominantly non-local part. This is evidenced by the strong stabilization of the 
ferromagnetic state for Mn and Cr pairs in GaAs at ~ 8 A separation. This interaction 
induces a spin-polarization of the host-like states. The direction of the spin-polarization 
depends on the relative energy position of the cation vacancy generated (host-like) states 
with respect to the impurity states. Furthermore, the band-theoretic description of Cr in 
GaP shows a partially occupied mid-gap band, and the wavefunctions associated with this 
mid-gap state are localized. Yet, even in the absence of free carriers, our total energy 
calculations predict a ferromagnetic ground state to be strongly stabilized, while no long 
range magnetic order is expected. 



II. EARLIER ELECTRONIC STRUCTURE CALCULATIONS 



There have been considerable earlier first-principles 



work on 



the electronic structure of these systems. One of the most well-studied systems is Mn in 



GaAs, which is found to be half-metallic 
in it has a net magnetic moment of 4 /z^ 
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121 1. The GaAs cell with one Mn atom 
, with part of the moment residing 



on the As neighbors of the Mn atom. The hole state is found to be strongly hybridized 
with the transition metal state, and has been referred to as a hybridized band of holes 
The Mn atom and its four nearest neighbors are found to account for most of the density 
of states at the valence band edge pj]. Since only the first shell of As atoms surrounding 
Mn are affected by the spin polarization of the Mn atom, the interactions are believed to be 



short-ranged 







The first principle results have been interpreted by Ref. as suggestive 
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of a d^/cf' electron configuration on Mn [6|. LDA+U calculations j3| have been used to 
obtain a description of the electronic structure consistent with photoemission. 

First-principle calculations have been used earlier to elucidate the magnetism in these 
systems. Mahadevan and Zunger developed a simple model of interaction of the cation- 
vacancy generated states with the transition metal states to understand how ferromagnetism 
results when Mn is doped into the chalcopyrite semiconductor CdGeP2. Sato and Katayama- 
Yoshida ^ have calculated the energy difference between ferromagnetic and the random 
alloy to determine which impurity could give rise to ferromagnetism. They found that at 
low concentrations V, Cr, and Mn doping in III-V stabilized the ferromagnetic state, while, 
Fe, Co and Ni doping stabilized the magnetically disordered state. Mirbt, Sanyal and Mohn 
|l2j | showed that the interaction of the transition metal impurity with the As dangling bond 
states could result in a spin polarization of the hybridized dangling bond states. The partial 
occupancy of these spin-polarized levels results in ferromagnetism. Sanvito, Ordejon and Hill 
B found an decrease of the spin-splitting of the valence band maximum of the GaAs host 
with impurity concentration. This is in contrast with what is expected from the mean-field 
Kondo Hamiltonian traditionally used to describe these systems Q]- They attribute the 
deviation to a breakdown of the mean-field approximation, while they say that the Kondo 
Hamiltonian isgood enough to provide a description of the ferromagnetism. Schilfgaarde 
and Mryasov lr\ have used the total energies obtained from first-principle calculations for 
different materials to extract exchange interaction strengths. They find a decrease of the 
exchan ge i nteractions with concentration which prompts them to suggest that the picture 
jl^ I15I. llfi. 12, 1^ of carrier-mediated ferromagnetism is not valid for these systems. 



We build on the current understanding of the electronic and magnetic properties that 
exists in the literature. However, we focus our calculations specifically on the examination 
of the features (i)-(iii) assumed as "input" to most model Hamiltonian theories. 

III. METHOD OF CALCULATIONS 

We have carried out first-principle electronic structure calculations using de nsity func- 
tional theory, within the momentum space total energy pseudopotential method j2J], using 
ultra-soft pseudopotentials [25] as implemented in the VASP j26[ code. The Ga pseudopo- 
tentials that we used for GaAs and GaP did not include the Ga 3d states in the valence. 
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While this is usually a good approximation for GaAs and GaP, it has been found that for 
GaN not retaining Ga 3d states in the valence leads to erroneous results for some physi- 
cal properties such as optimized lattice constant, cohesive energy We therefore used 
ultra-soft pseudopotentials which included Ga d states in the basis for GaN. We studied 
transition metal impurities - V, Cr, Mn, Fe, Co and Ni in 64 atom supercells of zincblende 
GaSb, GaAs, GaP and GaN. In Table I we compare the calculated lattice constants of pure 



III-V using PW91 GGA exchange functional 



28| with the experimental values 



29|. We have 



fixed the equilibrium lattice constant of the supercells at the calculated values of the pure 
host given in Table I. The basis sets had a cutoff energy for plane waves equal to 13.3 Ry 
for GaSb, GaAs and GaP and 29.4 Ry for GaN. We used a Monkhorst Pack grid of 4x4x4 
k-points which includes F. The cell-internal positions of the atoms were allowed to relax 
to minimize the forces. The equilibrium transition metal-to-As bond lengths in GaAs were 
2.47, 2.47, 2.48, 2.44, 2.36 and 2.34 A for V, Cr, Mn, Fe, Co and Ni respectively. 

The d partial density of states as well as the local moment at the transition metal were 
calculated within a sphere of radius of 1.2 A about the atoms and have been broadened with 
a gaussian of 0.2 eV full width at half maximum. The total energy differences between the 
ferromagnetic and the antiferromagnetic spin arrangements were computed for TM pairs at 
first and fourth neighbor separations for a parallel and an anti-parallel arrangement of their 
spins to determine whether a specific transition metal impurity resulted in a ferromagnetic 
state or not. 

and GGA 28j exchange ^ functionals, 
we consider the case of Co impurity in GaAs, where earlier LDA work 12] suggests a 
nonmagnetic ground state. Using the experimental lattice constant of 5.65 A for GaAs, we 
find that the GGA calculations lead to a magnetic ground state with a moment of 2 ^ib- 
The energy of this state is strongly stabilized by ~ 150 meV compared to the nonmagnetic 
state. Using LDA exchange functional we find that while the nonmagnetic state is stabilized 

n 

for a 2x2x2 Monkhorst Pack grid as used in the earlier work |12|, the magnetic state with 
the moment of 2 /x^ is stabilized by ~ 40 meV for a 4x4x4 Monkhorst Pack k-points grid. 
These observations are consistent with the fact that GGA calculations have a greater ability 
to stabilize a magnetic ground state than LDA calculation. For other impurities - Cr and 
Mn in GaAs, the LDA and GGA results are found to give the same ground state. We use 
the GGA exchange functional throughout this work. 
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The introduction of various transition metal impurities lead to defect levels in the band 
gap of the semiconducting host. We compute the formation energies of the transition metal 

or a defect comprising of atoms 



impurities in various charge states q. The formation energy 
a in the charge state q was computed using the expression 

AH^^\ef,fi) = E{a) - E{0) + J2rio.fil + q{E. + ep), (1) 

a 

where E{a) and E{0) are, respectively the total energies of a supercell with and without the 
defect a. rta denotes the number of atoms of defect a transferred in or out of the reservoir, 
while /i^ denotes their chemical potentials. 

Total energies: The total energies of the charged supercells were computed by compen- 
sating any additional charge on the impurity atom by a uniform jellium background and 
have been corrected for interactions between charges in neighboring cells using the Makov 
and Payne correction 3^. We use the static dielectric constant values - 15.69 for GaSb, 
12.4 for GaAs, 11.11 for GaP and 10.4 for GaN |33]. The quadrupole moment of the iso- 
lated defects was calculated as the difference between the moments of the supercell with the 
charged defect and that with the neutral defect. 

Transition energies: The defect transition energy e(g, q') is the value of the Fermi energy 
ep at which AH"''^(ef)=AH°''^' (ef). The zero of the Fermi energy is chosen as the valence 
band maximum E^ of the pure host. 

Chemical potential limits: As the reservoir supplying the atoms could be elemental solids, 
or compounds formed from the elements, we express /i^ as the sum of the energy of the 
element in its most stable structure /i^, and an additional energy /Iq, z.e /i^ = /i* + fia- The 
stable structures we considered for the elements were nonmagnetic body-centred cubic (bcc) 
for V, antiferromagnetic bcc for Cr, the antiferromagnetic face-centred cubic (fee) for Mn, 
ferromagnetic bcc for Fe, ferromagnetic hexagonal for Co, ferromagnetic fee for Ni and the 
nonmagnetic base-centred orthorhombic structure for Ga. 

The required ranges of //q are determined by fica < 0; ^itm < 0; fisb,As,p,N < (no 
precipitation of solid elements) and by the formation energies of the semiconducting host 
and competing binary phases formed between the elements of the semiconductor and the 
transition metal impurity. This could be the most stable NiAs phase of MnAs in the case 
of GaAs:Mn and the MnP phase of CrAs for GaAs:Cr. 

The energies E(q;), E(0), and /Xq, entering Eq.(l) are calculated within the density func- 



tional formalism discussed earlier. No correction for the band gap underestimation was 
made. Changing the k-point mesh from 2x2x2 to 4x4x4 changed the formation energies by 
~ 20 meV. We used a plane wave cutoff of 13.3 Ry for the calculations involving Mn in 
GaAs. Increasing the cutoff to 29.4 Ry, changed the formation energies by ~ 10 meV. 

IV. RESULTS OF DENSITY FUNCTIONAL CALCULATIONS 

We now divide the main features of the first principle calculations into three main entities 
introduced in Section I: 

A. The nature of the impurity- induced level in the gap 

B. The impurity-induced valence-band resonances 

C. The perturbed host VBM 

Then, in section V, we will provide a simple model that explains all of our numerical 
results qualitatively. 

A. The Nature of the impurity-induced level in the gap 

Figure 1 shows the transition-metal local DOS for V, Cr, Mn, Fe, Co and Ni in GaAs, 
projected into irreducible representations t2 and e and spin directions + and -. The VBM 
is at the zero of the energy. The GGA band gap of pure GaAs is found to be 0.3 eV; all the 
impurities V-Ni introduce levels into this band gap. We first discuss the nature of these gap 
levels, and then the circumstances how and when a hole is present in them. 

From Fig. 1 we see that the sequence of levels occupied for Cr in GaAs are t+, e+, t_ and t+ 
in the order of increasing energy. For a free atom one would expect levels of one spin channel 
to be filled up before levels of the other spin channel; the deviation that one observes here 
reflects solid state effects. The two sets of t+ and t_ levels that we find for each impurity are 
suggestive of bonding/ antibonding combinations arising from hybridization. We therefore 
determined the atoms on which each of the ^2 states are localized by computing atom- 
projected DOS. Bonding states with a large wavefunction amplitude on the TM site are 
referred to as "Crystal Field Resonances" (CFR) [J, whereas antibonding ^2 states with 
low contribution on the TM which are localized instead on the four nearest As atoms are 
referred to as the 'Dangling Bond Hybrid' (DBH). The full explanation of the genesis of 
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FIG. 1: The TM projected density of states for substitutional V-Ni impurities in GaAs evaluated 
in a sphere of radius 1.2 A for spin-up t+ (shaded region), spin down t_ (solid black line) and e 
(dashed line) symmetries. The zero of energy represents the valence band maximum of the host. 
The number of k-points used is 64. 

these states will be provided in Sec. V. We see that, 

1. Symmetry of gap levels and lowest unoccupied levels: Substituting Cr, Mn and Co 
in GaAs introduce levels with up-spin character and t2 symmetry in the band gap. These 
levels are partially occupied by 1, 2 and 2 electrons for neutral Cr, Mn and Co, respectively: 
CrO (tX), MnO (t^) and Co^ (tl). The levels introduced by V° (e^) and Fe° (tlel) are fully 
occupied. The first unoccupied levels have t+ and e_ symmetry for V and Fe respectively. 

2. d character of gap levels: The transition metal projected partial density of states for 
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FIG. 2: The accumulated charge Q within spheres of radius Ri and Ri+ARi about the TM 

impurity for Mn in GaN (filled squares), Cr in GaAs (open inverted triangles) and Mn in GaAs 
(filled circles) compared with the result for an electron gas (solid line), a is the lattice constant of 
the host supercell. 

different transition metal impurities in GaAs given in Fig. 1 indicates that the gap level/first 
unoccupied level is strongly d-\ike for the early transition metal impurities V and Cr, while 
for the heavier 3d elements e.g. Mn these levels have less d character. An increased d 
character of the gap level would imply increased spatial localization of the wavefunction in 
the vicinity of the impurity. 

3. Degree of localization of gap levels: We quantify the degree of localization by plotting 
in Fig. 2 the charge Q(Ri)= J^^^^ ijj'^r'^dr enclosed between concentric spheres with radius 
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Ri and Ri + A centered about the impurity atom. The integrated charge between the 
spheres is plotted as a function of Ri. For comparison, we show also the result expected for 
a homogeneous charge distribution (electron gas), where the charge density at any point in 
the cell is given by reciprocal volume y. We see that Q(Ri) for the TM impurities has little 
similarity to the results for an electron gas. Changing the impurity from Mn to Cr in GaAs, 
we see an increase in the charge density localized in the vicinity of the impurity. We find 
that till a radius which includes second neighbors of the TM atom, the integrated charge 
for Cr is higher than for Mn. Further, we find that the enclosed charge in the vicinity of the 
impurity atom is higher in GaN:Mn than in GaAs:Mn and the decay of the wavefunction is 
faster. 

4. The negative exchange splitting of the gap levels: Having established the identity of 
the gap levels, we now investigate their spin splittings. In Fig. 3 we have plotted the spin 
splittings of the CFR and DBH levels at F point for the impurities V-Co in GaAs obtained 
from an analysis of their eigenvalues/eigenfunctions. For V and Cr the spin splitting of the 
DBH levels is positive i.e. t^^^ states are at lower energies compared to t^^^ . However, for 
Mn, Fe and Co the splitting is negative with the t^^^ states at lower energies compared to 
t9^^ . A similar negative exchange splitting was observed earlier for the Te states in MnTe 
3^ . the Ce states in CeFe2 3^, and the Mo states in Sr2FeMo06 36]. The explanation 
was that the p states of Te in CdTe, the d states of Ce in CeFe2 and the d states of Mo in 
Sr2FeMo06 are sandwitched in between the 3d states of the transition metal atom. The p-d 
hybridization results in an exchange splitting of these states opposite in direction to that 
on the transition metal atom. Indeed we see from Fig. 1 that for the cases where the DBH 
states are bracketed by the spin split CFR states, the spin splitting is negative. 

5. Enhanced exchange splitting for t'-'^^ states: From Fig. 3, we see that the exchange 
splitting of the t'"^^ states is larger than that of the e'"^^ states for Mn, Fe and Co impurities. 

Having discussed the existence of impurity-induced levels in the band gap of the host 
semiconductor, we next discuss the location of these levels. 

6. Acceptor transitions for gap levels: Single particle LDA or GGA levels do not have 
any rigorous meaning. We thus calculate transition energies, e{q,q') which correspond to the 
value of the Fermi energy ep aX which the defect changes from a charge state q to q'. Table 
II provides the calculated and measured acceptor/donor transition energies for various 
transition metal impurities in GaAs. The calculated acceptor levels for Mn and Cr in GaSb, 
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FIG. 3: The Exchange-splitting at F point for t^^^, e^^^ (upper panel) and t^^^ (lower panel) 
states for V-Co impurities in GaAs. 

GaAs, GaP and GaN are plotted in Fig. 4 where the host band edges are ahgned according to 



their calculated unstrained valence band offsets [37[ . We see that as the electronegativity of 
the host crystal increases in the sequence GaSb GaAs GaP GaN, its bulk ionization 
energy (= position of VBM with respect to vacuum) increases. The acceptor level is thus 
farther away from the VBM of GaN than it is from the VBM of GaAs. Thus, GaN:Mn and 
GaP:Cr have more localized hole states whereas GaSb:Mn has more delocalized holes. This 
behavior, whereby the acceptor energy level does not follow the host valence band energy 
(as in the case of hydrogenic impurities) characterizes localized states 
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GaN 



GaP 



GaAs 



GaSb 



FIG. 4: The (0/-) acceptor transition energies for Cr and Mn impurities in GaN, GaP, GaAs and 
GaSb. The band edges of the host semiconductors are aUgned according to the LDA-calculated 
unstrained band offsets [s?!, and the gaps are the experimental values. 

7. Multiplet states and violation of isovalency rule: We use the level occupancies (Fig. 1) 
as well as the net magnetic moments that we obtain for different transition metal impurities 
in GaAs (Table III) to obtain the multiplet configuration describing the ground state. These 
are given in Table III. We also provide the multiplet configuration observed from experiment 

and find that there is agreement in all cases. It is interesting to compare the ground 
state multiplets of two isoelectronic cases - ZnSe:Fe^^ and GaAs:Co'^^ in their neutral charge 
states. In both cases we expect an electron configuration of d^. Normally, one would expect 
to find equal multiplets for isoelectronic cases jl| (isovalency rule). This expectation is based 
on the fact that we are looking at a low transition metal impurity concentration regime where 
basic crystal field theory ideas are expected to be sufficient to explain the observed ordering 
of energy levels. However, we find in GaAs:Co^+ the configuration '^T2, i.e the hole is in the 
t^^^ level (tS^^^ 

^'cFR+ ^CFR- ^DBH~ ^DBH+)^ whereas in ZnSe:Fe^+ we find ^E, i.e the 
hole is in the ccfr level (tc'FR+ ^cfr+ ^dbh- ^\>bh+ ^cfr-)- The reason for the difference 
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is that the stronger p-d hybridization for GaAs:Co pushes the t^^^ levels to higher energies, 
so that the e'^^^ levels are occupied first and the hole resides in the t^^^ level. We thus 
conclude that the isovalency rule is not applicable, and one cannot assume that the hole is 
in a "generic" d state. 

8. FM vs. AFM ground state and their relation to the symmetry of the gap levels: Having 
summarized the nature of the level induced in the gap by the introduction of the transition 
metal impurity, we now analyze when a ferromagnetic state is favored. In Table 111 we 
provide the energy difference between the ferromagnetic and antiferromagnetic energies for 
two TM atoms at nearest AEtvat and fourth neighbor AE47VAr fee positions in a 64 atom 
supercell of GaAs. We find that (a) when the level in the gap is fully occupied as in V° and 
Fe°, the favored ground state is antiferromagnetic. (b) When the level in the gap is partially 
occupied and has t2 symmetry as in Cr°, Mn", the ferromagnetic state is lower in energy. This 
is also the case for electron doped V~ in GaAs which is strongly ferromagnetic. Although 
Co*^ (^T2) also has a hole in the t2 level, the system is at the brink of a ferromagnetic-to- 
nonmagnetic transition, (c) When the level in the gap has e symmetry, as in the case of 
electron doped Fe~, the stability of the FM state is weaker. Evidently the symmetry of 
the hole carrying state strongly determines the magnetic order. We conclude that FM is 
stabilized strongly only when the hole resides in the level with t2 symmetry. (Note , viz. Sec. 
V, that in symmetry t2 states are strongly bonded to their neighboring atoms, whereas 
the lobes of the e orbitals point in between the nearest-neighbor atoms.) 



B. The impurity-induced valence-band resonances 

Figure 1 shows that in addition to the gap levels, the introduction of a transition metal 
atom gives rise to resonance levels that lie deep within the valence band of the host semi- 



one usually ignores the 



conductor. In most model Hamiltonian theories \U 
orbital degree of freedom of the transition metal impurity, and the presence of the impurity 
is included only as a localized spin of value 5/2. In our calculations we find that the de- 
gree of localization of such deep resonances (thus, the possibility of depicting them as local 
point-like spin) varies sharply with the position of the impurity in the periodic table. For 
heavier TM such as Fe and Mn, the deeper resonance level has significant TM d character 
(being Crystal-Field Resonances), while for the early TM impurities in GaAs, one finds that 
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the deeper t2 levels have significantly less TM character (being Dangling Bond Hybrids). 
This is discussed next. 

1. Anticrossing of the two t2 levels in different host materials: Level anticrossing is evident 
when keeping the impurity atom fixed, and, changing the host semiconductor. Considering 
the example of Mn, we find that by changing the host from GaSb to GaN, the DBH and 
the CFR exhibit anticrossing. This is not the only difference: We find that the exchange 
splitting of the DBH levels is in the same direction as the CFR levels (positive) in GaN:Mn, 
in contrast to GaAs:Mn. Further, in GaN:Mn the t'^^^ levels lie above the e'^^^ levels, 
unlike the case in GaAs:Mn. The reason is evident from Fig. 4 which shows that the VBM 
of GaN is much deeper than the VBM of GaAs. Since the free Mn^"*" ion has its d orbitals 
above the GaN VBM, but below the VBM of GaSB or GaAs, an anticrossing occurs along 
the GaN GaP GaAs GaSb series. This is illustrated in Fig. 5 which shows that in 
GaN:Mn for the up spin channel, the upper t2 is more localized than the lower t2, whereas 
in GaAs:Mn the localization sequence is reversed. This clarifies a confusion that existed in 
the literature regarding the question of whether the gap level is localized or not. Our 
result shows that the answer depends on the host. These results also clarify the nature of 
the acceptor transition for Mn in different materials. GaN:Mn can be viewed as a d^-like 
case since its configuration is e^(jpj^_^t^pp^j^ (hole in t'^^^), and the (0/-) acceptor transition 
is from a Mn configuration d^ to d^. On the other hand for Mn in GaP, GaAs and GaSb 
we have the configuration {d^ + hole), i.e. e^Qp^^j^ tcFR+ ^%bh- ^'dbh+ the acceptor 
transition is from a Mn configuration of (d^+hole) to d^. 

2. Occupancy of the valence band resonances and comparison with photoemission: Table 
HI gives the calculated occupancy of the crystal field resonances of the impurities in 
GaAs (in square brackets). These levels are found to have a configuration of "c?^" for Cr, 
Mn and Fe. Experimentally the position of these levels can be detected by valence band 



photoemission j^l • By suitably tuning the photon energy so that the photo- ionization cross- 
section is maximum for the TM-related states, an electron can be ionized from these deep 
CFR level. Koba.asM . Q used .e.nant valence band p«oen...on and ..owed 
that the the CFR levels for Mn in GaAs are located at - 4 eV. A direct comparison of 
the position of these levels with the single-particle density of states calculated for GaAs:Mn 
places these energies at - 2-3 eV. The LDA error in the position of these states is because 
of the self-interaction correction (SIC) that places these energies too high [4^. As pointed 
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Energy (eV) 



FIG. 5: The up and down spin Mn d projected partial density of states evaluated within a sphere 
of radius 1.2 A for a Mn impurity in GaN (upper panel) and GaAs (lower panel). The number of 
k-points used is 64. 

out earlier j^] for the 3d states in II-VI's, the experimental result should be compared with 
the total energy difference between the configurations d'^ and d^ and not with the bare single 
particle eigenvalues. Alternatively, the LDA error can be empirically corrected by using 
the simplified LDA+U version of SIC. In Fig. 6 we plot the Mn t^^^ partial density of 
states as a function of U for GaAs:Mn. As U increases, the position of the Mn t'"^^ level 
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FIG. 6: Mn projected partial density of states evaluated within a sphere of radius 1.2 A for 
on-site Coulomb interaction strengths U=0, 6 and 10 and 15 for Mn in GaAs. The number of 
k-points used is 4x4x4. The inset shows the variation in E(FM)-E(AFM) for two Mn atoms at 
nearest neighbor positions. 

is pushed deeper into the GaAs valence band. Agreement with XPS for the t'^^^ being at 
E^-A cV occurs for U ~ 2 eV. At the same time, the Mn character of the DBH state at ep 
(not shown) decreases with increasing U. As the DBH hole becomes more delocalized, the 
Efm-^afm stabilization energy (insert to Fig. 6) is reduced; a "host-hke hole" obtained for 
unphysically large U leads to nearly vanishing FM stabilization energy. Clearly, the picture 
of "host like hole" is invalid for GaAsrMn, since for the U that leads to agreement with XPS 
the DBH hole is still localized to some extent, whereas for very large U, when the hole is 
delocalized, there is no ferromagnetism. 
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GaP:Cr GaPrCr GaP 

(a) up-spin (b) down-spin (c) host 




FIG. 7: The band dispersions in (a) up and (b) down spin channel for GaP:Cr compared with (c) 
the host supercell. The thickness of the hnes represents the Cr weight in the bands. 

C. The perturbed host VBM 

Having studied the impurity-induced levels in the gap and deep in the host valence 
band, we next examine the perturbation of the host states, especially the host valence band 
maximum by the presence of the impurity atom. Figure 7 shows the up and down spin 
band dispersions for 3 % Cr doped GaP supercell (panels (a) and (b)). The band dispersion 
of the GaP host without the impurity has been provided in panel (c) for comparison. The 
thickness of the lines depicting these bands has been made proportional to the Cr d character 
of the states. We see that Cr introduces a new band within the band gap of GaP. In a band- 
theoretic picture, this system is metallic, with the Fermi energy within the impurity band. 
Interestingly, (1) the host band dispersions are significantly altered by the presence of the 
impurity. In particular the VBM is found to have significant TM d character for the 3% Cr 
concentration represented by the supercell. (2) A Cr-induced spin-sphtting of the valence 
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FIG. 8: Up (solid line) and down (dashed line) spin projected partial density of states for As atoms 
labelled 1-4 in Mn substituted GaAs evaluated within spheres of radius 1.2 A using 64 k-points. 
The positions of the As atoms (filled circles) with respect to the Mn atoms is shown in the inset. 

band maximum is observed. Effects (1) and (2) suggest that the host VBM is sufficiently 
perturbed by the transition metal. 

Another way of detecting perturbations in the host bands is to examine the host projected 
DOS of the system containing the impurity. In Fig. 8 we have plotted the As p partial density 
of states projected onto different As atoms labelled 1-4 for a GaAs supercell containing 2 
Mn atoms. The As atom labelled 1 has one Mn nearest neighbor, while the As atom labelled 
2 has two Mn nearest neighbors. The As atoms show a strong polarization which increases 
with the number of Mn neighbors. The As atoms labelled 3 and 4 which are far away from 
the Mn atoms show a reduced polarization. 

To pictorially see the perturbation in the VBM states, we compare in Fig. 9 the wave- 
function squared along two chains in the (110) plane for the valence band maximum of the 
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pure GaAs host (panel (a)) as well as the VBM {i.e. state below DBH) of the system with 
the Mn impurity in the up (panel (b)) and down (panel (c)) spin channels. The upper chain 
in panels (b) and (c) contains the perturbing Mn impurity. The perturbation of the VBM 
in the presence of the impurity can be assessed by comparing the perturbed charge density 
for each spin channel with the unperturbed charge density of the host lattice. We find that 
the perturbations are significant in the chain containing the Mn atom, while in adjoining 
chains, the perturbation is limited in extent. Further the perturbations are stronger in the 
up spin channel than in the down spin channel. 

To evaluate the 3c?-induced spin-splitting in the VBM, we reference the up and down 
spin VBM eigenvalues of the impure system to the corresponding VBM of the pure host 
semiconductor. This is done by aligning the average potentials on Ga atoms far away from 
the impurity for the two systems. The presence of the impurity band with t2 symmetry 
above the VBM for Mn and Cr impurities complicates the identification of the valence 
band maximum. We associate the highest occupied triply degenerate eigenvalues at F point 
with the impurity band, and the next deeper set as E^. The shift with respect to the 
pure host A£;+=£^+(GaAs:Mn)-£;,„ (GaAs) and AE-=E-{GaAs:Mii)-E., (GaAs) is given 
in Table IV for the impurities V, Cr and Mn in different host semiconductors. We find that 
the perturbation of the host VBM is smaller in the down spin channel compared to the 
up spin channel. This is consistent with what we find from the charge density plotted in 
Fig. 9. The spin splitting of the valence band maximum depends strongly on the transition 
metal impurity and host semiconductor combination. Considering the case of impurities in 
GaAs, we find that while the spin splitting associated with the introduction of V is only 
0.06 eV, it increases to 0.39 cV for Mn. Keeping the impurity fixed (Mn), and varying the 
semiconductor host (GaAs to GaN), we find the splitting decreases from 0.39 eV to 0.1 eV. 
The small valence band splittings in the case of V in GaAs as well as Mn in GaN compared 
with that for Mn in GaAs is because of the larger energy separation between the interacting 
t2 states in the former cases compared to the latter. 

D. Summary of the electronic structure as obtained by density functional 

Wc arc now in a position to examine whether the physical picture of the electronic 
structure of 3d impurities as assumed in model Hamiltonian theories (reviewed in Sec. I) is 
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Wavefunction at Valence band maximum at r 




FIG. 9: The wavefunction squared at the valence band maximum for (a) the pure host GaAs 
(b) for GaAs:Mn in the up spin channel and the (c) down spin channels. The lowest contour 
corresponds to 0.0015 e/A^ and each contour is 1.6 times larger. 
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consistent with first-principles calculations (outlined in Sec. IV). 

(i) The nature of the TM-induced hole state: A 3d impurity in a III-V semiconductor 
generates two sets of states with t2 symmetry, and one set of states with e symmetry in each 
spin channel. While one set of t2 states are localized on the TM atom (CFR), the other are 
localized on the host anion atoms next to the impurity (DBH). These states CFR and DBH 
exhibit anticrossing for a fixed TM as a function of the host anion GaN — GaP — ^ GaAs — >■ 
GaSb, or for a fixed host as a function of the impurity V —> Mn. The localization of the hole 
state decreases as we move from Mn in GaN to Mn in GaP, and then to Mn in GaSb. Not 
all impurities introduce holes. In GaAs, V" and Fe" have no hole; Cr", Mn° and V~ have t2 
holes; and Fe~ has an e hole. In all cases, however, the hole is non-hydrogenic, manifesting 
significant admixture of 3d character and showing deep acceptor levels whose energies do not 
follow the host VBM. This implies that the neglect of the short-range part of the impurity 
potential and the consequent expansion of the acceptor wavefunction in terms of a single 
host wavefunction are questionable. The effective mass of the hole state is therefore different 
from that of the host, as observed in recent experiments [4^. The exchange splitting of the 
CFR states is different for the t2 states from that for the e states. While the splitting for 
the e states is larger than that for the t2 states for V and Cr in GaAs, the order is reversed 
for Mn, Fe and Co. This reversal in the order of the spin splitting of the CFR states is 
accompanied by a reversal in the sign of the spin splitting of the DBH states. The identity 
of the hole state - both the symmetry as well as the character depends on the impurity-host 
combination. While the hole carrying orbital for Fe in ZnSe has e symmetry, the hole is 
found to be located in an orbital with t2 symmetry for the isovalent doping of Co in GaAs. 

(ii) The nature of the host VBM: The introduction of the transition metal perturbs the 
valence band of the host crystal. We find this perturbation to be large when the state in 
the gap has t^^^ character. This is because the effective coupling is larger since the DBH 
states have strong host character. We find that the VBM is spin-split in the presence of 3d 
impurity, and that the VBM in the up-spin channel is perturbed more strongly than in the 
spin down channel in the presence of the impurity. 

(iii) Ferromagnetism and symmetry: Impurities with fully occupied DBH-like ^2 gap states 
such as V°, Fe*^ show antiferromagnetism. Partial occupation of the t^^^ as in Cr°, Mn° 
or V~ show ferromagnetism. Partially occupied e-like level like in Fe~ show weak or no 
ferromagnetism. 
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(iv) Ferromagnetism and hole localization: Using LDA+U as an artificial device to explore 
the consequences of delocalized host-like-hole states we find (insert to Fig. 8) that in this 
limit there is reduced ferromagnetism. 

We find that despite the well known GGA-LDA band gap error, as well as the under- 
estimation of the location of deep CFR states due to SIC, these first principle calculations 
provide us with the correct spin multiplets. LDA+U changes some details (CFR locations), 
but does not alter the basic picture emerging from GGA/LDA when the hole is DBH-like 
as in GaAs:Mn. 

V. SIMPLE MODEL OF THE ELECTRONIC STRUCTURE OF 3d IMPURITIES 
IN GAAS 

A. The Model 

Most of the results of the density functional study of the electronic structure of 3d im- 
purities in III-V's (Sec III) can be captured by a simple model. The electronic structure of 
substitutional 3d in III-V semiconductors can be understood as arising from the interaction 
of the host cation vacancy (= anion dangling bonds) with the crystal-field and exchange 
split orbitals of a 3d ion. 

(a) The dangling bonds for a column III cation vacancy Vni'- A cation vacancy for a 
column III element gives rise to a fully-occupied s-like ai level located deep in the host 
valence band, and a partially-occupied p-like ^2 level located just above the host valence 
mnd maximum, with wavefunction amplitude localized primarily on the neighboring atoms 
4^ . This is evident from the wavefunction squared of the Ga vacancy dangling bond state 
shown in the (110) plane in Fig. 10. The neutral vacancy V^^j has a deficiency of 3 electrons, 
i.e the orbital configuration is a{t\{p), where p denotes its major orbital character. Spin- 
polarization splits this t2{p) vacancy level into spin-up [t+(p)] and spin-down [t-(p)] states, 
but the splitting is small (90 meV at the F point) on account of the rather delocalized nature 
of these pure host dangling bond orbitals. 

(b) The crystal-field split TM 3d orbitals: The tetrahedral crystal-field of the zincblende 
host splits the TM d levels into e{d) and t2{d), with e below t2 in the point-ion limit the 
crystal-field (CF) splitting of the ion is denoted by AcF{'t2-G)- Exchange interactions further 
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Ga-vacancy dangling bond states in GaAs at r 




(b) dn spin (unoccupied) 



FIG. 10: The wavefunction squared for the (a) up and (b) down spin Ga-vacancy generated 
danghng bond states with t2 symmetry in pure GaAs. The lowest contour corresponds to 0.0017 
e/A^ and each contour is 1.6 times larger. 
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FIG. 11: The schematic energy level diagram for the levels (central panel) generated from the 
interaction between the crystal-field and exchange-split split levels on the 3d transition metal ion 
(left panel) with the anion dangling bond levels (right panel), when the TM d levels are energetically 
deeper than the dangling bond levels. 

split these levels into spin-up (+) and spin-down (-), with exchange-splittings Ax{e) =[e_(d)- 
e+(d)] and A,{t) =[t_(d)-t+(d)]. 

The energy levels of a cation-substituted TM in a III-V semiconductor can be thought of 
as the result of hybridization between the anion dangling bonds generated by a column III 
cation vacancy Vjjj [(i) above], and the crystal-field and exchange split d levels of a TM ion 
placed at the vacant site [(ii) above]. There are two limiting cases: When the 3d levels are 
well below the host cation dangling bonds (e.g. Mn in GaAs, Fig. 11), or when the 3d levels 
are well above the host cation dangling bonds (e.g. V in GaAs, Fig. 12). The dangling bond 
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FIG. 12: The schematic energy level diagram for the levels (central panel) generated from the 
interaction between the crystal-field and exchange-split split levels on the 3d transition metal ion 
(left panel) with the anion dangling bond levels (right panel), when the TM d levels are energetically 
shallower than the dangling bond levels. 

states are shown on the right hand side of Fig. 11 and 12, while the crystal field and exchange 
split TM d levels are shown on the left hand side of Fig. 11 and 12. The levels generated after 
hybridization are shown in the central panel. The t2{p) levels of the anion dangling bond 
hybridize with the t2{d) levels of the transition metal. In contrast, the e{d) level of the TM 
ion remains largely unperturbed since the host does not have localized e states in this energy 
range, available for significant coupling. Considering the examples of Mn in GaAs and GaN, 
we plot the charge density of the e and the t2 states in the (110) plane in Fig. 13. It is 
evident that the e states for Mn in GaAs are essentially nonbonding, while in GaN, as a result 
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of the reduced Mn-Mn separation, there is weak interaction between the Mn atoms. The 
hybridization in the i2-channel creates bonding, transition-metal locahzed "Crystal-Field 
Resonances" (CFR): t'^^^, t^^^, as well as the host-anion locahzed antibonding "Dangling 
Bond Hybrids" (DBH): t^^^ , t^^^ , whereas the e-channel creates the non-bonding e'j^^^ 
and e^^^ states. This model explains the existence of two sets of and t_ levels that we 
found in Fig. 1. The available electrons for occupation of these levels are N=(n-l)-|-6 for a 
d"s^ transition metal atom (three electrons are used to complete the anion dangling bond 
state to leaving d"~^ at the transition metal ion.) For GaAs:V (Fig. 12) the ordering 

of levels after hybridization is 

^DBH ^ ^DBH ^ ^CFR ^ ^CFR ^ ^CFR ^ ^CFR 

with increasing energy. Hence for V we have N=(n-l)+6=8 electrons occupy the DBH and 
CFR levels. Thus, V has the configuration {t%BH+^DBH-^CFR+)^ ^ ^^^^ ^^S- 1 
Table HI. 

The order of levels for 

GaAs:Mn is {tcpji^ < e^FR+ < '^dbh- < ^'dbh+)j 

GaAs:Fe is (tc'Fii+ < ^cfr+ < ^dbh- < 

GaAsrCo is {tcpji_^ < e^FR+ < ^cfr- < ^dbh- < t'DBH+) 

(Fig. 11). The number of electrons (n-l)+6 is 10, 11 and 12 for Mn, Fe and Co respectively. 
This agrees with Fig. 1 showing that Mn and Fe in GaAs have the ordering of levels shown 
in Fig. 11, with fully-filled t^^^ and e%^^ levels and 2, 1 and holes in the t^^" level. 
By an analysis of the density of states obtained within our first principle calculations, we 
have determined (Table HI) the energy minimizing orbital configurations for the transition 
metal impurities V, Cr, Mn, Fe and Co in GaAs in fully relaxed configurations. The first 
unoccupied orbital for each impurity has been indicated in boldface in Table HI. The simple 
model of Figs. 11, 12 gives the same result. 



B. Qualitative consequences of the simple model 



1. Level anticrossing: The model explains how the hopping interaction between the 
t2 states on the transition metal impurity with the cation-vacancy states generates a pair 
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Mn-induced state in GaAs and GaN at r point 
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FIG. 13: The wavefunction squared of Mn induced (a) e*^^^ in GaAs:Mn, (b) e*^^^ in zinc blende 
GaN:Mn, (c) t^^^ in GaAs:Mn, (d) t^^^ in zinc blende GaN:Mn, (e) t^^" in GaAs:Mn and (f) 



iDBH 



in zinc blende GaN:Mn. The lowest contour corresponds to 0.015 e/A^ and each contour is 



1.6 times larger. 



of t2 states in each spin channel. The bonding-antibonding character of these states is 
determined by the relative separation of the interacting levels as well as their interaction 
strengths. Hence as depicted in Fig. 14, one could by a suitable choice of the TM impurity 
change the character of the gap levels. When the orbital energy of the 3d ion hes below the 
host dangling bond, we have a "CFR-below-DBH" situation, illustrated in Fig. 11. In this 
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Changing 3d impurity in GaAs 
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GaAs:l\/ln,Fe GaAs:V, Cr 

FIG. 14: The schematic plot of band anticrossing between the two t2-Hke levels in GaAs for 
different 3d impurities. 

case one has CFR states in the valence band of the semiconductor while the gap levels are 
more delocalized with dominant weight in the dangling bonds. This is the case for GaAs:Fe, 
Mn, and Co. Conversely, when the orbital energy of the 3d ion lies above the host dangling 
bond, we have the " CFR-above-DBH" situation, illustrated in Fig. 12. In this case the gap 
level is CFR-hke. 

While Fig. 1 illustrates anticrossing when changing the 3d atom, but keeping the host 
fixed, e.g. GaAs, Fig. 5 suggests that there is also anticrossing when keeping the 3d 
atom fixed e.g. Mn, but changing the host crystal GaN — ^ GaP GaSb. Indeed, as the 
anion dangling bond becomes deeper (GaSb GaAs —>■ GaP GaN) or the TM d level 
becomes shallower (Mn — V), there will be a DBH-CFR anticrossing, resulting in the level 
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ordering shown in Fig. 12. The CFR-above-DBH is exemphfied by GaAs:V (see Fig. 1) 
which has a shallow d level and by GaN:V, Cr and Mn which have a deep semiconductor 
VBM js^. For GaAs:V, this level anti-crossing results in a reduced /S.x{t2) exchange splitting 
compared to the cases where this crossing did not occur (Cr and beyond). For Mn in III-V's 
we see a reversal between hole-in-DBH case for GaAs:Mn with \t\e\t^_e^_]cFR {^Xt^-SjoBH 
configuration to the case of hole-in-CFR [t^e^t° e° ]cfk {t%-tt)DBH for GaN:Mn. As the 
CFR is much more localized than DBH, this is refiected by a change in acceptor level depth. 

2. Negative exchange splitting in the DBH manifold and low-spin configurations: For 
Co,Fe and Mn in GaAs the special position of the Ga- vacancy level t=p(p) between the 
exchange-split TM ion levels t+(d) and t-(d) (Fig. 11), results in a hybridization- induced 
exchange splitting of the DBH states (Fig. 3), opposite in direction to the splitting on the 
CFR levels (Fig. 3). This results from the fact that t^^^ is pushed down by t^{d)-t-{p) 
coupling more than the t+{d)-t^{p) coupling^ pushes t^^^ up. Such a ne ^ive exchange 
splitting was also observed in bulk MnTe CeFa ^ and SraFeMoOg As seen m 
Table III, the negative exchange splitting leads to the DBH orbital configuration {t^t^) and 
(t+t"^) for Mn and Fe, respectively. This corresponds to a low-spin configurations for Mn 
with a moment of 4. Co, on the other hand, has a configuration of [t^e^e'^]cFR (t^'tX)DBH, 
so fi=2fiB- Note that the moment //=4 for GaAs:Mn is a consequence of A^, < 0: the e'^t^ 
CFR levels give a spin of 5/2, and the t:^'^^-below-t:^'^^ gives a configuration t^t"^, so the 
total spin is 4/2, and /i=4. Had t^^^ been above t^^^ , we would have had the configuration 
t+t^ with S=3, and /i=6. As a result of the negative exchange splitting of the DBH states, 
one finds that the exchange splitting for the more delocalized t'"^^ levels is larger than that 
for the more localized e'"^^ levels. 

3. Migration of the d-holes into the DBH acceptor states of GaAs:Cr, Mn and Co: The 
configuration of the trivalent TM ion corresponds to d'^ for Cr and d"^ for Mn in GaAs 

and thus to 2 and 1 hole respectively, relative to d^. However, since for these impurities the 
t^^^ and e^^^ levels are deeper in energy than the t^^^ levels (Fig. 11), it is energetically 
favorable for the d holes to "fioat" into the higher t^^^ levels. As a result, we find that 
GaAs:Cr (d^) has a configuration [ti]_e^]ci?/j DBH with two holes in the DBH, while for 
GaAs:Mn (c/^) we find [t^e'^]cFR{t^t'+)DBH with one hole, and GaAs:Fe (d^) has a closed- 
shell configuration of [t^e'^e'^]cFR(t+t^)DBH- One expects Co (d^) to have a configuration 
[t^e'^el_]cFR(t-t^)DBH, i-G. a ground state multiplet with a hole in e*"^'^. However, 
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we find that the lowest energy configuration for Co is [t\e\e'^_]cFR{t^t'^)DBH with a hole 
in the DBH level i.e. ground state multiplet ^T2. The migration of the holes from the 
TM-locahzed deep CFR levels into the dangling-bond manifold for Cr, Mn and Co in GaAs 
creates partially occupied p-d hybridized states at the Fermi level with significant delocalized 
As p character. 

4. Perturbation of the valence band maximum of the host: The valence band maximum 
of the pure host has t2 symmetry. Hence, they can interact with the states with the same 
symmetry on the dangling bonds as well as the TM ion. The consequent spin sphtting 
of the VBM must depend on the relative separation of the levels involved as well as the 
coupling strength. Our analysis for different TM-host combinations suggests that when 
the gap level has t^^^ character the spin-splitting of the VBM is large (~ 0.4 eV for 
Mn in GaAs). Changing the host semiconductor from GaAs to GaN, increases the energy 
separation between the dangling bond levels generated by the Ga vacancy and the valence 
band levels. Consequently the perturbation of the valence band levels resulting in the 
observed spin-splitting is smaller (~ 0.1 eV for Mn in GaN). This VBM spin-sphtting of 
the impure system has been traditionally used to estimate the exchange interaction strength 
Jpd between the hole and the transition metal atom. We see that the spin splittings of the 
VBM of the impure system is grossly underestimates the spin-splitting of the impurity band 
(DBH,CFR) as seen in Figs. 1,3 and 5. 

VI. CONCLUSIONS 

We have analyzed the basic electronic structure of 3d transition metal impurities in III- 
V semiconductors. Wc find that the introduction of a 3d impurity is accompanied by the 
introduction of a pair of states with ^2 symmetry in addition to nonbonding states with e 
symmetry. Not all 3d impurities introduce holes. The basic symmetry and character of the 
hole state depends strongly on the semiconductor-impurity combination. We find that the 
hole has significant 3d character. We have constructed a microscopic model which captures 
the basic aspects of the electronic structure of transition metal impurities in semiconductors. 
The elements of this model are the relative separation of the dangling bond and transition 
metal levels, the p-d hybridization strength and the crystal-ficid and exchange splittings of 
the transition metal levels. We model a change in these interaction strengths by changing 
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the semiconducting host, keeping the transition metal impurity fixed, - Mn. We find that 
while the hole introduced by Mn has significant 3d character in GaN, it is more delocahzed in 
GaAs. The symmetry (e vs. ^2), the character (DBH vs. CFR) as well as the occupancy of 
the gap level, determine the magnetic ground state favored by the transition metal impurity. 
When the hole has dominantly host character, an exchange splitting is induced on the hole 
states which is opposite in direction to that on that the transition metal atom. The nature 
of the exchange coupling Jpd that exists between between the transition metal atom and 
the hole comes out automatically from such a microscopic model. The perturbation of the 
host valence band is not directly related to the couphng strength Jpd. When the hole has 
primarily DBH character, one finds the perturbation of the host valence band is larger. The 
basic picture that emerges from our first-principles calculations could be used to replace the 
more naive model-Hamiltonian treatments which have assumed a host-like-hole picture, an 
unperturbed valence band and that the spin of the hole couples to the spin of the TM via 
a local exchange interaction. 

This work was supported by the U.S. DOE, Office of Science, BES-DMS under contract 
no. DE-AC36-99-G010337. 
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TABLE I: Comparison of GGA optimised lattice constants with experiment for the pure host 



System 


experiment a (in A) 


GGA PW 91 a(in A) 


GaN 


a=4.49; 4.51 


a=4.53 


GaP 


a=5.45 


a=5.489 


GaAs 


a=5.65 


a=5.728 


GaSb 


a=6.10 


a=6.18 
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TABLE II: Impurity formation energies with and without Makov-Payne charge corrections for 
the acceptor transitions for 3d impurities in GaAs. Experimental transitions have been given in 
brackets for comparison. 



System 


Formation energies 
(no charge correction) 


Formation energies 
(with charge correction) 


GaAs:V (q=0) 
GaAs:V (q=+l) 
GaAs:V (q=-l) 
GaAs:V (q=-2) 


1.22 +IJ,Ga - fJ'V 
1.32 +IJ,Ga - IJLV + 
2.05 +HGa - tiV - 
3.18 +IJ,Ga - fJ'V - 2eF 
(0/-) = 0.83 eV 
(-/2-) = 1.13 eV 


1.22 
1.41 
2.15 
3.56 

(0/-) = 0.93 eV 
(-/2-) = 1.41 eV 


GaAs:Cr (q=0) 
GaAs:Cr (q=+l) 
GaAs:Cr (q=-l) 
GaAs:Cr (q=-2) 


1.61 +I^Ga - f^Cr 
1.47 +HGa - IJ-Cr +^F 
2.02 +HGa - IJ-Cr -CF 

2.81 +/LtGo - A*Cr -2ei7 
(0/-) = 0.41 (0.74) cV 
(-/2-) = 0.79 (1.57) eV 


1.61 
1.56 
2.115 
3.20 

(0/-) = 0.51 (0.74) eV 
(-/2-) =1.09 (1.57) eV 


GaAs:Mn (q=0) 
GaAs:Mn (q=+l) 
GaAs:Mn (q=-l) 


1.04 +HGa - fJ'Mn 
1.15 +HGa - IJ'Mn +^F 
1.13 +IJ,Ga - IJ-Mn -^F 

(0/-) = 0.09 (0.11) eV 


1.04 
1.24 
1.23 

(0/-) = 0.19 (0.11) eV 


GaAs:Fe (q=0) 
GaAs:Fe (q=+l) 
GaAs:Fe (q=-l) 
GaAs:Fe (q=-2) 


1.79 +HGa - fJ'Fe 
1.83 +/XGo - fJ'Fe +eF 
2.10 +HGa - fJ'Fe -^F 
2.82 +fJ,Ga - fJ-Fe '^^F 

(0/-) = 0.31 eV 
(-/2-) = 0.72 eV 


1.79 
1.92 
2.21 
3.27 

(0/-) = 0.42 eV 
(-/2-) = 1.06 eV 


GaAs:Co (q=0) 
GaAs:Co (q=+l) 
GaAs:Co (q=-l) 


1.84 +IJ.Ga - IJ-Co 
1.90 +HGa - IJ-Co +^F 
1.92 +^Ga - ^^Co -CF 

(0/-) = 0.08 (0.16) eV 


1.84 
1.99 
2.01 

(0/-) = 0.17 (0.16) eV 
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System 




Formation energies 


Formation energies 






(no charge correction) 


(with charge correction) 


GaAs:Ni (q= 


=0) 


1.73 +/iGa - fJ'Ni 


1.73 


GaAs:Ni (q= 


=+1) 


1.76 +HGa - fJ'Ni +eF 


1.86 


GaAs:Ni (q= 


-1) 


1.86 +HGa - fJ'Ni -CF 


1.96 


GaAs:Ni (q= 


=-2) 


2.27 +fiGa - fJ-m -2ei? 


2.68 






(0/-) = 0.13 (0.22) eV 


(0/-) = 0.23 (0.22) eV 






(-/2-) = 0.41 (1.13) eV 


(-/2-) = 0.72 (1.13) eV 



TABLE III: The calculated energy-minimizing configuration for neutral substitutional 3d impuri- 
ties in GaAs. CFR states are given in square brackets and DBH states in round brackets. Boldface 
letters denote the first unoccupied orbital. Also shown are ground state multiplet and in parenthe- 
ses, local moment ^loc within a sphere of radius 1.2 A for isolated impurities. The last 2 columns 
give the total energy difference AE between FM and AFM spin arrangements of TM pairs at first 
(NN) and fourth neighbor (4NN). Asterisk denotes the configuration with lowest energy. 



TM 


Configuration 


Multiplet 




AE4NN 






{fJ-loc) 


(in meV) 


(in meV) 


Ni 




(0.53) 


+2.85 


+4.3* 


Co 


[tlelt'Lel] {tHD 


(1.58) 


-9.6* 


-22.6 


Fe 


[^e^to eO] {tHD 


Ml (3.27) 


+298* 


+205 


Mn 


[tleltH°_] {tHD 


(3.75) 


-247* 


-227 


Cr 


[tleltH^-] {tHD 


^Ti (2.99) 


-315* 


-258 


V 


[tO e2.t0 eO ] {tHD 


=^^2 (1.84) 


-40 


+31* 


V- 


[tO.Se^tOgO] (t3^3) 




-204 




Fe- 


[tDltH^j^] {tHD 




+259* 
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TABLE IV: Spin-up and spin down eigenvalues at F for 3d impurities in GaN, GaP and GaAs, 
referenced to the valence band maximum of the pure host. A corresponds to the spin splitting 
between up and down spin states. 



System 


E+ (eV) 


K (eV) 


A (eV) 


GaAs:V 


-0.14 


-0.08 


0.06 


GaAs:Cr 


-0.26 


-0.08 


0.18 


GaAsiMn 


-0.47 


-0.08 


0.39 


GaP:Cr 


-0.22 


-0.08 


0.14 


GaP:Mn 


-0.39 


-0.09 


0.30 


GaN:Cr 


-0.07 


-0.03 


0.04 


GaN:Mn 


-0.13 


-0.03 


0.10 
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